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ABSTRACT 

As an extension of the method by Bahcall et al. (2004) to investigate the time dependence 
of the fine structure constant, we describe an approach based on new observations of forbidden 
line multiplets from different ionic species. We obtain optical spectra of fine structure transitions 
in [Ne III], [Nc V], [0 III], [01], and [SII] multiplets from a sample of 14 Seyfert 1.5 galaxies in 
the low-z range 0.035 < z < 0.281. Each source and each multiplet is independently analyzed to 
ascertain possible errors. Averaging over our sample, we obtain a conservative value a 2 (i)/a 2 (0) 
= 1.0030±0.0014. However, our sample is limited in size and our fitting technique simplistic 
as we primarily intend to illustrate the scope and strengths of emission line studies of the time 
variation of the fine structure constant. The approach can be further extended and generalized 
to a " many-multiplet emission line method" analogous in principle to the corresponding method 
using absorption lines. With that aim, we note that the theoretical limits on emission line ratios of 
selected ions are precisely known, and provide well constrained selection criteria. We also discuss 
several other forbidden and allowed lines that may constitute the basis for a more rigorous study 
using high-resolution instruments on the next generation of 8 m class telescopes. 



Subject headings: galaxies: active - Atomic data 

1. Introduction 

The variation of the fine structure constant a 
=e 2 /Aireohc is of fundamental interest in cosmol- 
ogy. However, if there is a variation of the fun- 
damental constants by time, the effect will be 
very small. Recent laboratory measurements us- 
ing m Yb+ by Pcik ct al. (2004) give an up- 
per limit of 2xl0~ 15 yr~ 1 resulting in a change 
of a of the order of 10 -5 in 10 Gyrs. In or- 
der to measure this effect using astronomical ob- 
servations, they have to be performed with high- 
rcdshift objects. Therefore the lines to be studied 
need to manifest themselves as strong features in 
the spectrum. The most common methods are 
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measurements of high-redshift Lyct forest metal 
absorption lines. Among the first ones was the 
"alkali-doublet method" (e.g. Bahcall et al. 1967) 
with transitions from the singlet ground level to 
the two fine structure doublets in alkali-like sys- 
tems, such as the recent study of C IV and Si IV 
systems using a QSO spectra from UVES (Mar- 
tinez Fiorenzano et al. 2003). The most exten- 
sive body of work in recent years is described by 
Murphy et al. (2003) (and references therein), 
who have considerably extended absorption line 
studies to a relatively large number of multiplets 
in heavy atomic systems, including the iron group 
elements. Their "many multiplet method" repre- 
sents 'an order of magnitude increase in precision' 
over the alkali doublet method and uses Lymana 
forests of Quasars (Dzuba et a. 1999a, b; Webb et 
al. 1999, 2002; Murphy et al. 2001a,b, e.g.). 

Recently, as an alternative Bahcall et al. 
(2004) described an extensive analysis of forbid- 
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den fine structure lines of the well known nebular 
[O III] doublet at 5006.84 and 4958.91 A. These 
lines have the great advantage that they are ex- 
tremely bright and nearly ubiquitous in the optical 
spectra of H II regions in most sources. The ap- 
proach by Bahcall et al. (2004) involved the study 
of archived spectra of 3,814 quasars from the large 
database of the Sloan Digital Sky Survey (SDSS 
York et al. 2000). They devised elaborate al- 
gorithms based on stringent selection criteria to 
search for and obtain a standard spectral sample 
of 42 quasars, as well as alternate samples based 
on variants of the selection criteria. Among these 
criteria are fits to the line profiles and line ra- 
tios. Since both lines originate within the same 
upper level, their line profiles must be similar and 
their line ratios depend only on intrinsic atomic 
properties, the energy differences between the fine 
structure levels and corresponding spontaneous 
decay Einstein A-coefficients, i.e. independent of 
external physical conditions such as the density, 
temperature, and velocities. Bahcall et al. (2004) 
also noted that other similar forbidden multiplcts 
of [Ne V] and [Ne III] could be exploited for emis- 
sion line studies, but they arc likely to be weaker 
than the [O III] by an order an magnitude in 
typical quasar spectra. 

The quasar absorption line many- mult iplet 
method has clear advantages: because it mea- 
sures the Lya forest of high-rcdshift systems at 
observed optical wavelengths its look-back time 
into the Universe's history is long. It also uses 
many line pairs to measure any time dependence 
of a and provides good statistics to minimize sys- 
tematic errors of the measurements of an indi- 
vidual line pair. However, the method has some 
disadvantages: It is observationally expensive, re- 
quiring very high resolution, like A/AA >50000 
using high-resolution Echcllc spectrographs and 
long exposure times on large telescopes (e.g. Mur- 
phy et al. (2003)). Also the line pairs measured 
by this method are from different atomic stages 
suggesting that they are not necessarily formed 
in the same regions. Another disadvantage is 
that the large number of absorption lines in Lya 
forest systems make it difficult for the correct 
identification of the lines; many lines are blended 
requiring multi-Gaussian fits to the line blends 
in order to determine the wavelengths. Based 
on high-resolution Keck spectra of several differ- 



ent samples of quasars, Murphy et al. (2003) 
report a 'highly significant' value of Aa/ao = 
-0.57 ± 0.10 x 10~ 5 . In contrast the statistically 
invariant result by Bahcall et al. (2004) yields a 
value of 0.7 ± 1.4 x 10~ 4 . An extensive discussion 
of relative problems and advantages of various 
methods is given, among others, by Bahcall et al. 
(2004) and Levshakov (2003). 

The optical [O III] forbidden emission lines 
AA4959, 5007A studied by Bahcall et al. (2004) 
are relatively easy to measure in most AGN, (ex- 
cept in Narrow-Line Seyfert 1 galaxies (NLSls) 
in which these lines can be very weak and con- 
taminated by strong Fell emission (e.g. Boroson 
& Green 1992; Grupe 2004; Sigut & Pradhan 
2003, 2004). The observations can be performed 
with lower resolution Spectrographs and smaller 
telescopes. Bahcall et al. (2004) presented a sam- 
ple of 44 AGN carefully selected from the SDSS 
and measured the wavelength shifts between the 
[OIII]AA4959,5007A lines. The disadvantage of 
this method is that it only uses one line pair 
and it can only be used in the optical wavelength 
range for objects with redshifts z<0.8 before the 
[OIII]A5007A line is shifted out of the observable 
optical wavelength range and NIR spectroscopy is 
required for objects of higher rcdshift. 

In this paper we generalize the emission line 
method outlined by Bahcall et al. (2004) by in- 
corporating, in principle, the main advantage of 
the many-multiplct absorption line method. How- 
ever, there are other key differences with Bahcall 
et al. (2004). We carry out new observations of 
quasars and Seyfert 1 galaxies, rather than analyze 
existing data. This obviates the need for complex 
search routines that forms the bulk of the analysis 
of the SDSS data by Bahcall et al. (2004) . We ex- 
pect better accuracy since objects are pre-selected 
from previous observations for good quality opti- 
cal spectra to enable optimum analysis. Since the 
redshifts are known a priori, we know the approx- 
imate positions of the well known forbidden lines. 
Theoretically known line ratios (discussed in the 
next section) offer an additional indication of the 
quality of the spectra. 

This many- mult iplet emission line method not 
only allows us to measure more than one line pair 
per object and gives a more secure result, it also, 
once established, opens the possibility of observing 
objects of higher rcdshift in the optical wavelength 
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range. Another purpose of our observing run was 
to test whether it is feasible to use a 2m class tele- 
scope with a medium resolution spectrograph to 
get enough accuracy to perform a measurement of 
the time variation of the fine-structure constant a. 
While the availability of observing time at large 8- 
10m class telescope is very limited, many institu- 
tions have access to 2-3m small/medium size tele- 
scopes. In this work we focus on laying out the 
framework for general emission line analysis, po- 
tentially leading to future studies using the many- 
multiplct emission line method, with a relatively 
small sample. However, we discuss a variety of 
further extensions. 

This paper is organized as follows: in § 2 we de- 
scribe the theoretical background of the fine struc- 
ture separation of the LS multiplets, in § 3 we 
describe the spectroscopic observations at MDM 
observatory, in § 4 the results are presented, and 
discussed in § 5 usually in the context of systemat- 
ical errors in the wavelength calibration. Finally, 
we note a few main features of our method in the 
concluding section 6. 

2. Theory 

Relativistic electron interactions lead to fine 
structure separation of LS-coupling multiplets. 
The level energies may be expanded in terms of 
a non-relativistic term and relativistic terms in 
powers of a (Bethc & Salpeter 1977). Although 
the exact formulation is predicated on the Dirac 
theory (Grant 1996), Breit formulated the gen- 
eralized interaction for non-hydrogenic systems 
that leads to the so called Breit-Pauli approxima- 
tion for the total electronic Hamiltonian as a sum 
of one-body and two-body operators, referred to 
as the non-fine structure and fine structure terms 
(Drake 1996; Frose-Fisher 1996). The leading rel- 
ativistic terms, such as the spin-orbit coupling, are 
of the order of a 2 ; higher order terms are orders 
of magnitude smaller. Forbidden lines arise from 
transitions between LS states of the same elec- 
tronic configuration of an atomic system, usually 
the ground configuration. Fine structure compo- 
nents of a forbidden LS multiplet are separated 
by AE ex a 2 . The time dependence of a there- 
fore should manifest itself in different energy or 
wavelength separations at different cosmological 
epochs. 



2.1. Forbidden Lines and a (t) 

Bahcall et al. (2004) show that, for two lines 
from a forbidden multiplet, the ratio 



R(t) 



X 2 — Ai 



(1) 



A 2 + Ai ' 

at a cosmological time t is related to the ratio at 
the present epoch t = by 



R(0) a 2 (0)' 



(2) 



where the RHS is a measure of the variation 
in a as a function of the 'cosmological look-back 
time' t. Wavelength separation between a pair 
of well defined forbidden lines may thus be used 
as a 'chronometer' provided its value at an earlier 
cosmological epoch differs from that at present. 

As mentioned, Bahcall ct al. (2004) had al- 
ready suggested the possibility of extending this 
study with other emission line pairs, such as 
[Ne V]AA3346,3426A and [Ne III]AA3869,3968A. 
The left panel of Figure 1 shows the schematic di- 
agram of four ionic species ( [Nelll] , [NeV] , [OIII] , 
and [OI]) and the two lines 1 D 2 ^ 3 P2 an d 
1 D 2 Pi- It might be noted that while [O I] 
and the isoelectronic ions C-like [O III] and [Ne V] 
have the ground level 2p2( 3 P ), the O-like [Ne III] 
has a different atomic structure with the ground 
level 2pi( 3 P 2 ). In addition we consider the well 
known nebular doublet lines [Oil] and [SII] due to 



transitions D° 



CO 

°3/2- 



The atomic 
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transitions within the ground configuration cor- 
respond to magnetic dipole (Ml) and the much 
smaller electric quadrapole (E2) interactions. The 
wavelengths and atomic transitions of the [Nelll] , 
[NeV], [OIII], [OI], and [SII] lines are given in 
Table 1. 

2.2. Line Ratios 

Subsequent analysis of different emission line 
multiplets rests on the basic property of the ratio 
of lines originating in the same upper level. For 
such a 3-level system the theoretical emissivity ra- 
tio of line intensities is 



LR = 



N 3 A 31 hv 31 



(3) 



N 3 A 32 hv 32 

The upper level population is the same for both 
lines, and the R t depends only the intrinsic atomic 
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parameters, the transition rates or A-values, and 
the energy differences. Generally, for forbidden 
lines, the A-values arc obtained from sophisticated 
theoretical calculations, whereas the energies are 
available from laboratory measurements. While 
the prime resource for both quantities is the on- 
line National Institute of Standards and Technol- 
ogy (NIST) compilation (www.nist.gov), some of 
the NIST data for A-values is not quite up-to-date. 
The differences may be slight but quite significant. 
Tabic 2 shows the line ratios for the 5 ions and re- 
spective multiplets using both the NIST data and 
the most recent calculations. 

In their work Bahcall et al. (2004) quote the 
line ratio for the [O III] AA4959,5007A line using 
NIST values as 2.92, as opposed to their measured 
value of 2.99 ± 0.02. However, the 2.92 is the ra- 
tio of the A-values alone, not the cmissivity ra- 
tio above; taking account of the energy differences 
the NIST tabulations yield 2.89. While this is in 
even worse agreement with the measured ratio, the 
most recent calculations of A-values by Storey & 
Zeippen (2000) yields a significantly better the- 
oretical value of 2.98, in agreement with Bahcall 
et al. (2004). The crucial difference between the 
Storey & Zeippen (2000) results and previous cal- 
culations (such as by Galavis et al. 1997) is taking 
account of higher order rclativistic corrections to 
the magnetic dipole Ml operator (Drake 1971; 
Eissner & Zeippen 1981). The above discussion 
emphasizes the need for extremely accurate atomic 
calculations including all relevant relativistic ef- 
fects, as well as a well optimized configuration in- 
teraction expansion, both of which are important 
to obtain precise A-values for fine structure tran- 
sitions. 

The actual value of the forbidden line ratio may 
be of decisive importance in spectral identification 
and analysis of observational datasets, and they 
indicate the extent of line blending and signal- 
to-noise ratio. It is therefore instructive to ex- 
amine further (with a view towards future work 
as well) the line ratios of interest in this work. 
For the forbidden lines in the left panel of Fig- 
ure 1 the dominant contribution is from the Ml 
operator (the E2 contribution is about 3 orders 
of magnitude smaller). In the LS coupling limit, 
as the magnetic interaction goes to zero, the ratio 
of the line strengths S( 1 D 2 - 3 P 2 )/S( 1 D 2 - 3 Pi) 
= 3 (Storey & Zeippen 2000). The A-values 



however involves the energy differences as well, 
i.e. A( 1 D 2 - 3 P 2 )/A( 1 D 2 - 3 P0 = 3 E{ 1 D 2 - 3 
P 2 )3/E{ 1 D 2 - 3 Pi)3. The line ratio (Eq. 3) de- 
viates from the value of 3, in accordance with the 
magnitude of the magnetic interaction. It is highly 
fortuitous that for O III this ratio is in fact 3, as 
noted by Bahcall et al. (2004) and (Storey & Zeip- 
pen 2000). But for other ions there is significant 
deviation, as seen in Table 2. We do not give line 
ratios for the [SII] lines because their line ratios 
depends strongly in the gas temperature (e.g. Os- 
terbrock 1989). We derive this ratio for all three 
ions using 4 different datasets available in litera- 
ture. While all sets of line ratios agree to within 
few percent, we envisage that this level of differ- 
ence could turn out to be crucial in analyzing large 
datasets, with line ratios measured and calculated 
to better than 1% accuracy. 

In addition to the ions with ground state fine 
structure levels discussed above, there are well 
known pairs of lines in ions with excited state 
fine structure levels, such as the doublet forbidden 
lines [O II] and [S II] at AA(A) = (3728.8/3726.0) 
and (6716.5/6730.8) respectively (right panel of 
Figure 1), observed from many H II regions and 
AGN (Ostcrbrock 1989; Pradhan 1976). While 
collisional coupling leads to variations in line ra- 
tios dependent on electron density, the limiting 
values of these ratios are precisely known. For ex- 
ample, the low density line ratio LR (N e — > 0) 
= 1.5 for both [O II] and [S II]: the ratio of the 
statistical weights of the upper levels in associated 
pair of transitions ( 2 D° /21 2 D° /2 — > 4 S° /2 ). The 
high density limit is given by 

_ g( 2 D° /2 ) A(( 2 D° /2 - 4 S° /2 ) 
LR{Ne °° ) " 9?BZ /2 ) A(^ /2 -^ /2 )< 

(4) 

and is = 0.35 for [O II] (Zeippen 1987) and 0.44 
for [S II] (Mendoza and Zeippen 1982). Studies 
of ionized gaseous nebulae show no known cases 
of deviations from these 'canonical' limits (Wang 
et. al. 2004). Whereas the wavelength separation 
between the [O II] doublet is possibly too small to 
be resolved in AGN, the [S II] doublet is easily ob- 
served and resolved, as in the present exploratory 
work. 
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3. Observations and data reduction 

We observed a small sample of 14 Seyfert 1.5 
galaxies with the 2.4m Hiltncr telescope at MDM 
observatory at Kitt Peak, Arizona, for 5 nights 
from 2003-10-16 to 2003-10-21 to cover the [NeV], 
[Nelll] and [OIII] lines — the blue part of the op- 
tical spectrum — and 3 nights starting 2004-10-10 
to observe the [OI] and [SII] doublets towards the 
red part of the spectrum. (Table 3). All spec- 
tra were taken with the OSU CCDS spectrograph 
with the 600 grooves/mm grating in first order. 
The slit width of the 2003 run was 1 correspond- 
ing to a spectral resolution in FWHM of about 2A. 
The weather conditions during the whole run were 
excellent with mostly photometric conditions. The 
slit position was normally at E-W for the spectra 
in the 4000-7000A range, but all spectra in the 
blue were observed in N-S direction to compen- 
sate for refraction losses in the earth's atmosphere. 
During the 2004 run the weather conditions were 
clear but suffered from rather bad seeing. Most 
observations during the 2004 run were performed 
with slit widths of 1.5 and 2.0 . The slit was ori- 
ented in N-S direction for all 2004 observations. 

The total exposure times per spectrum arc 
given in Table 3. Each spectrum consists of 2- 
4 single spectra which were combined after going 
through all data reduction steps. For each individ- 
ual spectrum a wavelength calibration spectrum 
was taken. For the wavelength calibration we used 
a Hg lamp for the 3400-4300A range, Xe for the 
4100-5100A range, Ne for the 5100-6000A range, 
and Ar for the 4900-5900A and all observation 
with A > 6000 A. For the flat field correction the 
CCDS spectrograph only allows internal flats. To 
perform a flat field correction we used an average 
of 10 flats. We used standard stars BD+28-4211, 
Feige 110, G191 B2B, and Hiltner 600 for the flux 
calibration. The data reduction was performed 
with the ESO MIDAS data reduction and analy- 
sis package version 01FEBpll.4. 

The wavelength measurements were performed 
by fitting a Gaussian at 80% of the line peak in or- 
der to avoid possible line asymmetries that can oc- 
cur at lower parts of an emission line. Line fluxes 
were measured by integrating over the whole line 
using the MIDAS command integrate/line. In or- 
der to estimate the error of the data reduction, the 
data were reduced and analyzed independently by 
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4. Results 

Figure 2 displays the optical spectra of the 
[NeVO, Nelll], and [OIII] line regions of the ob- 
jects sorted by RA as given in Table. 3. As shown 
in Figure 2 in most of the sources the lines are 
clearly present and have enough S/N that their 
wavelengths and fluxes can be measured. The 
spectra of the [OI] and [SII] lines taken in Oc- 
tober 2004 are shown in Figure 3. In most cases 
the [SII] lines arc clearly present. However, the 
[OI]A6363 line in most cases is too noisy to allow 
accurate measurements of the wavelength. 

Table 4 lists the line flux of the [OIII]A5007A line 
and the line ratios. In general, the [NeV]A3426 
and [NeIII]A3869 line are about l/10th of the 
[OIII]A5007 line. Table 5 lists the observed wave- 
lengths of the [NeV], [NcIII], [OIII], [OI], and [SII] 
lines. 

For the low-redshift AGN in our sample with 
a maximum look-back time of about 2 Gyrs we 
would expect a maximum change of a of 4xl0 -6 
regarding the the laboratory measurements of 
Peik et al. (2004). Therefore the expected ratio 
a 2 (t)/a 2 (t) maximum is 1.000008. Every devia- 
tion from this value gives us a handle on the qual- 
ity of our measurements. Table 6 lists the ratios 
of a 2 {t)/a 2 {Q) which are equal to the R(t)/R(0) 
ratio, with R(t) = (A2(t) - Al(t))/(Al(t) + A2(t) 
and R(0) are given by the laboratory wavelengths 
as R(0)([NeV])=1.1814, R(0)([NeIII])=12.597, 
R(0)([OIII])=4.80967, R(0)([OI])=5.011971, and 
R(0)([SII])=1.0690657 xlO" 3 (for all lines). Ta- 
ble 6 also summarizes the mean values of a 2 (t) /a 2 (0) 
for each object as well the sample averaged value. 
In general, we find results with the lowest error 
estimates from the [Nelll] lines owing to their 
wide wavelength separation, whereas the values 
obtained from the [SII] lines suffer in precision 
from their small splitting. Clearly, the lines with 
the highest SNR offer the best possibilities to con- 
strain their centers, making the original approach 
based upon the [OIII] doublet so effective. 

Even though we are not attempting to measure 
a redshift dependence of the fine structure con- 
stant a by redshift due to the low redshift of the 
AGN in our sample, Figure 4 displays the redshift 
vs. a 2 (t) / a 2 (0) diagram. As expected we do not 
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see any significant deviation from 1.0000. 

5. Discussion 

As mentioned before, it is beyond the scope of 
this paper to measure the time dependence of the 
fine structure constant a with adequate precision. 
The redshifts of the objects in our sample only 
cover a range between 0.034-0.281. The look-back 
time of an object of z=0.28 is of the order of ~ 3 
Gyrs. With the upper limit of a possible change 
of a measured by Peik et al. (2004) of Aa = 
2 x 10~ 15 yr _1 we would expect an upper limit of 
Aa = 6 x 10 -6 , which implies that a two orders of 
magnitude improvement is needed over what we 
can achieve from our current data set. 

For measuring any redshift dependence of the 
fine structure constant a , estimates of the ac- 
curacy in determining the wavelengths of the 
line doublets arc absolutely crucial. Evaluating 
a 2 (t)/a 2 (0) we do not need to rely on absolute 
wavelength calibrations, but relative wavelengths 
ratios. However, the issues that are important 
are: First, the stability of the wavelength calibra- 
tion between two lines in a line pair, which can 
be separated up to «100A in the case of [NcIII], 
and second, the ability to determine the line cen- 
troids. While the first issue is relatively safe for 
the [OIII], [OI] and [SII] lines, because their sep- 
aration is rather small and the wavelength cali- 
bration is very secure due to the large number of 
calibration lamp lines, the situation is more dif- 
ficult for the [Nelll] and [NeV] lines. Not only 
that the separation between the two lines in the 
[Nelll] and [NeV] line pairs are about 100A, also 
the wavelength calibration in the blue using the 
CCDS spectrographs suffers from a lack of enough 
calibration lamp lines. The CCDS spectrograph 
only has a Hg lamp providing just 5 calibration 
lines throughout the blue/UV wavelength range. 
Certainly, a Th/Ar lamp could improve our abil- 
ity to provide for a more stable calibration in this 
regime. The situation for these line pairs would 
also improve for objects of higher redshift as the 
observed features shift into the longer wavelength 
regime where the calibration is better. Further- 
more, we indicate that the abundant night-sky 
lines upwards from 6500 A could be used as addi- 
tional calibrators to measure small-scale flexures 
as already indicated by (Bahcall et al. 2004). 



Measuring the center of the lines depends on 
a variety of factors: the spectral resolution and 
dispersion of the instrument, the pixel sampling, 
the signal-to-noise ratio in the line and the line 
shape that often deviates substantially from a 
simple Gaussian. With the setup used for our ob- 
servations, the resolution is A/AA=2000 and the 
dispersion yields 0.79 A per pixel. A low S/N also 
introduces an error in the wavelength measure- 
ment, because noise changes the shape of the line 
and causes the estimate of the line center to shift. 
Clearly, lines with high SNR will provide the best 
results regarding this aspect. 

An important aspect of performing a high-precision 
analysis study is the sample selection. Bahcall et 
al. (2004) had to be very critical about the shape 
of the [OIII] lines of the sources in his sample 
reducing the number of good targets from about 
1000 to about 40. This was necessary for the 
SDSS sources which were reduced and analyzed 
automatically. For our small sample, we could 
work on each source manually and had hoped to 
obtain reasonably good lines from prior knowl- 
edge. Our sources were selected to be X-ray hard 
which tend to have stronger [OIII] emission than 
soft X-ray selected AGN (Grupe et al. 2004). For 
many of the sources optical spectra were avail- 
able from the ROSAT Bright Source Catalogue 
(www.aip.de/~aschopwe/rbscat/rbscat.html) which 
gave us some indication of the expected line 
shapes. Furthermore, by only fitting a Gaussian 
line to the narrow part of the emission lines, we 
tried to minimize the effect of asymmetric line 
shapes as much as possible. 

The identification and analysis of lines is greatly 
facilitated by the fact that the emission line multi- 
plcts chosen in our study have well constrained line 
ratios. Therefore observed line pairs with ratios 
outside the theoretical limits can be safely ruled 
out if they are blended or suffer from instrumental 
effects. However, it is essential that the relevant 
Einstein A-coefficients be accurately calculated, 
particularly taking account of all relativistic ef- 
fects. Whereas such data are available for the ions 
and transitions considered in this work, high preci- 
sion atomic calculations are needed before expend- 
ing the study to more complicated atomic species. 
If these calculations are carried out then emission 
line studies of other ions should be possible. 
Taking these factors into account individually for 
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each line pair, we have estimated an error bud- 
get for each line centroid measurement which is 
listed in table 5. On average, we thus believe to be 
able to constrain an individual line center to about 
0.6A, a value substantially higher than (Bahcall et 
al. 2004) who estimate their precision to 0.05 pix- 
els or about 0.06 A at 5000 A. The effect on the 
error budget of a(i) /a(0) - R(t)/R(0) depends 
on the actual wavelength splitting of the line pairs, 
and is most pronounced for the lines with small 
separations. We achieve the best values for the 
[NcIII] lines, but even in that case the error on 
an individual measurement of R(t) /R(0) is never 
below 1.0 x 10 -3 and usually of the order 7 x 10~ 3 . 

For the time being the primary aim of this 
study is to lay the groundwork for future studies 
based on the many emission-line method with 8- 
10 m telescopes to explore the high-z regime of 
faint objects. For example, the Large Binocular 
Telescope (LBT) is slated to have two spectro- 
graphs, the Multi-Object Double Spectrograph 
(MODS 2 Osmer et al. 2000) in the optical 0.3-1.0 
/^m range, and the other LUCIFER 3 in the J,H,K 
bands. Present observations were divided into two 
parts, one focusing on ions in the blue side ([NcV], 
[NcIII], and the other on ions in the red side of 
the optical spectrum ([OI] and [SII] ) , with [OIII] 
in the middle. As the highcr-z objects become 
accessible with, say, the LBT, pairs of lines from 
these ions would move into the range from MODS 
into the J,H,K bands covered by LUCIFER. This 
would enable a natural and logical extension of 
the present studies with LBT. 
With the predicted capabilities of LUCIFER, 
MODS at the LBT, we estimate the error bud- 
get of an individual line measurement of [NcIII] 
at a redshift of z~2.5 to be ~0.2 A which allows 
us to constrain a(t) 2 /a(0) 2 to 8 x 10~ 4 even with 
our much simpler fitting approach than (Bahcall 
et al. 2004) used for their sample. Implementing 
their detailed analysis for the determination of the 
line centroid and carefully choosing a reasonably 
sized sample with prior line-shape knowledge, we 
will be able to push the error limit for an averaged 
a(t) 2 /a(0) 2 to ~ 10~ 5 " 6 , comparable to the limit 
reached in recent absorption line studies and thus 
providing an interesting alternative. 

2 http: / /www. astronomy.ohio-state.edu/LBT/MODS/ 

3 http: / /www. lsw.uni-heidelberg.de/projects/Lucifer/indcx. html 



In summary, the emission line multi-multiplct 
method has certain advantages: it combines the 
multi-multiplet absorption line method, having 
many line pairs and getting good statistics of the 
measurements of a 2 (t)/a 2 (0), with the ease of 
identifying lines and measuring the wavelengths 
of the emission line pairs. It is observationally 
relatively inexpensive to get good spectra of the 
[OIII] and [Nelll] line pairs. However, our experi- 
ence with the current data set shows that by us- 
ing a 2m class telescope in order to archive better 
S/N of the [NcIII] and [NcV] lines more than the 
typical 1-2 hours of observing time that we spent 
have to be invested. Is this project still suitable 
or 2-3m class telescopes? In principle, yes. How- 
ever it becomes more and more challenging for ob- 
jects with higher redshifts which are fainter and 
therefore require much longer integration times. 
Because a possible time dependence of the fine 
structure constant a can only be measured from 
quasars with redshifts of z=3 or higher this type 
of high-precision measurements is limited to large 
telescopes with medium- or high-resolution NIR 
spectrographs only. However, because the reso- 
lution does not have to be as high as for the ab- 
sorption line method, the exposure time per object 
will be in the order of an hour only. Nevertheless 
smaller telescopes can build up the fundamentals 
which can be extended by larger telescopes. With 
a number of 8-10 m class telescopes available in 
the future it should be feasible to request alloca- 
tion for such a project. 

6. Conclusion 

We have demonstrated the application of a 
method for studying time variation of the fine 
structure constant based on the analysis of many 
emission line multiplets, with the following salient 
features: 

• Extension of the [O III] multiplct method 
by Bahcall et al. (2004) to a 'many multiplet 
method' of well known forbidden line multiplets 
in the optical rest frame. 

• Our best measurements archive errors in 
a 2 (i)/a 2 (0) in the order of 10 -3 . 

• In contrast to Bahcall et al. (2004) whose 
work based on a search of the SDSS database, the 
present work involved new observations of selected 
Seyfert 1 galaxies specifically targeted to obtain 
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several emission line multiplcts. Two separate ob- 
servation runs were made, focused on the blue and 
the red sides of the optical spectrum. 

• Best suited are sources with strong very nar- 
row NLR lines. However, this requirement be- 
comes more challenging for high-rcdshift AGN. 
High rcdshift AGN tend to have central black holes 
with masses in the order of 10 9 to lO lo M0(e.g., Di- 
etrich & Hamann 2004; Vestergaard 2004). Due 
to the well-known relation between the black hole 
mass and the bulge stellar velocity distribution, 
the Mbh-c relation (e.g., Ferrarese Merritt 2000; 
Gebhardt et al. 2000), high-redshift quasars will 
tend to have broader NLR emission lines than low- 
rcdshift AGN which tend to have smaller black 
hole masses. Counteracting this problem is the 
growing wavelength separation for increasing red- 
shifts. 

• If a Thorium calibration lamp is not available 
to observe at wavelength A <4100A, it is recom- 
mended to use only AGN with a redshift z>0.24 
to observe the [NeV] lines with appropriate wave- 
length calibration. 

• Analysis of observed pairs of lines showed 
the necessity of high accuracy theoretical atomic 
calculations in order to obtain line ratios which, 
in principle, can be ascertained a priori. 

• The next generation of 8- 10m class of tele- 
scopes should be able to achieve the required preci- 
sion of Aa(t) ~ 10 -5 " -6 to make a more definitive 
prediction. 

We would like to thank Axel Schwope (AIP) for 
making his ROSAT Bright Source Catalog avail- 
able to us. We would also like to thank Bob Barr 
and his crew at MDM observatory for their techni- 
cal support. This work was supported in part by a 
grant from the astronomy division of the National 
Science Foundation. 
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Table 1 

Rest-frame wavelength and atomic transitions of the line doublets 



Line 


A [A] 


Atomic transition 


[NeV] 


3345.86 


1 D 2 - 3 P 1 


[NeV] 


3425.86 


1 D 2 - 3 P 2 


[Nelll] 


3868.75 


^2 - 3 P 2 


[Nelll] 


3967.46 


X B 2 - 3 Pi 


[OIII] 


4958.91 


!D 2 - 3 Pi 


[OIII] 


5006.84 


^2 - 3 P 2 


[01] 


6300.304 


l B 2 - 3 P 2 


[01] 


6363.776 


r D 2 - 3 Pi 


[SII] 


6716.440 


2r\o 4 co 
^5/2 °3/2 


[SII] 


6730.816 


2r\o 4 co 
U 3/2 D 3/2 



Table 2 

Line Ratios For [O III], [Ne V] And [Ne III] Using Different Atomic Data. 



Ion 


AE( 1 D 2 - 3 P>) 


&E( l D 2 - 3 Pi) 


A( 1 D 2 - 3 P 2 ) 


A( 1 D 2 - 3 Pi) 


LR 


O III 


0.18195 


0.18371 


1.8105(-2) 


6.212 (-3) 


2.89 a 




?? 




1.96(-2) 


6.74 (-3) 


2.88 fe 




11 


11 


2.041(-2) 


6.995 (-3) 


2.89 c 




11 


11 


2.042(-2) 


6.785 (-3) 


2.98 d 


Ne V 


0.26592 


0.27228 


3.82(-l) 


1.38(-1) 


2.70 a 








3.65(-l) 


1.31(-1) 


2.72 b 




11 




3.499(-l) 


1.252(-1) 


2.73 c 




11 




3.501(-1) 


1.221(-1) 


2.80 d 


Ne III 


0.23548 


0.22962 


1.703(-1) 


5.24(-2) 


3.33 a 




11 




1.71(-1) 


5.42(-2) 


3.24 b 




11 




1.73(-1) 


5.344(-2) 


3.32 c 




11 




1.708(-1) 


5.413(-2) 


3.24 d 



Notes: a - NIST Compilation, b - Pradhan and Peng Compilation (1995), c - Galavis et al. (1997), d - 

Storey and Zeippen (2000). 
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Table 3 

Observation log, observing times are given in minutes 



# 


Object 


RA-2000 


DEC-2000 


B-mag 




z 


[NcV], [Nelll] 


T obH 
[OIII] 


[OI], [SII] 


1 


PG 0026+129 


00 29 14 


+ 13 16 04 


15.41 





14537 


60 


30 




O 




nn x,a 

UU Qt OZ 


1 OC oc on 


15 42 


n 
u 




fin an 


ou 


fill 
DU 


3 


RX J0334. 4-1513 


03 34 24 


-15 13 40 


15.43 





.03478 


90 


40 


1211 


4 


RX J0337. 0-0950 


03 37 03 


-09 50 02 


17.0 





.28074 


110 


(iO 


120 


5 


RX J0354. 1+0249 


03 54 09 


+02 49 30 


16.3 





03536 


80 


(ill 


120 


6 


RX J0751. 0+0320 


07 51 00 


+03 20 41 


15.2 





09914 


30 


25 


90 


7 


MS 0754+393 


07 58 00 


+39 20 49 


14.36 





09533 


60, 60 


20 


120 + 40 


8 


RX J0836. 9+4426 


08 36 59 


+44 26 02 


15.6 





25427 


90 


30 


75 


9 


MS 2128.3+0349 


21 30 53 


+04 02 30 


16.34 





08600 


120 


(iO 


120 


10 


PKS 2135-147 


21 37 45 


-14 32 55 


15.91 





.20048 


80, 60 


40 


90 


11 


RX J2256. 6+0525 


22 56 37 


+05 25 16 


16.2 





06529 


120 


40 


60 


12 


PG 2304+042 


23 07 03 


+04 32 57 


15.44 





04265 


00 


45 


80 


13 


RX J2325. 9+2153 


23 25 54 


+21 53 16 


15.9 





12033 


80 


(ill 


120 


14 


PG 2349-014 


23 51 56 


-01 09 13 


15.7 





17404 


60, 80 


50 


75 



Table 4 

Line fluxes (observer's frame). 



# 


Object 


F 1 
" [OIII] 4959 


p 

" [OIII] 5007 


p 

" [Nc V] 3346 


p 

" [NcV] 3426 


p 

^ [NcIII] 3869 


p 

-*^[NcIII]3967 


p 

F [OI] 6300 


p 

" [OI] 6364 


p 

" [SII] 6716 


p , 

" [SII] 6730 


1 


PG 0026 


304±30 


1014±40 


20± 10 


83± 20 


90± 20 


30± 10 





— 


— 


— 


2 


Pvj OOoz 


o2U± oU 


12yu± DU 


10± 5 


44 i 10 


lyO± 30 


45± 20 






16± 8 


25± 10 


3 


RXJ0334 


174± 30 


490± 30 






80± 5 


20± 10 


68± 8 


20± 10 


90± 15 


20± 5 


4 


RXJ0337 


126± 20 


400± 20 


10± 5 


50± 10 


64± 8 


28± 5 


45± 10 








5 


RXJ0354 


130± 6 


413± 15 


6± 4 


26± 10 


55± 10 


18± 4 


8.7± 2.0 


0.75± 0.7 


10. 9± 1.0 


8.4± 1.2 


(i 


RXJ0751 


540± 20 


1800± 100 


16± 5 


80± 10 


120 ±20 


50 ±30 


7 ±4 








7 


MS 0754 


980± 40 


3790± 100 


34± 10 


250± 50 


380± 30 


500 ± 30 






110± 30 


96± 4 


8 


RXJ0836 


390± 50 


1200± 60 


54± 10 


125± 10 


240± 20 


110± 20 


25± 10 


2± 2 






9 


MS 2128 


140± 15 


480± 20 


10± 5 


40± 10 


50± 10 


24± 4 


21± 12 


1.5± 1.0 






10 


PKS2135 


360± 10 


1190± 40 


9± 3 


50± 7 


55± 10 


20± 7 










11 


RXJ2256 


250± 10 


670± 30 




70± 10 


117± 20 


75± 50 


3.5± 1.5 




9± 2 


5± 3 


12 


PG 2304 


175± 25 


640 ± 30 




15± 10 


53± 10 




26± 6 


1.5 ±1.0 


7.5± 3.0 


11. 2± 4.0 


13 


RXJ2325 


146± 15 


520± 20 




24± 8 


42.5 ±5.0 


19± 7 


25± 5 








14 


PG 2349 


60± 20 


207± 30 






17. 5± 7.0 


3.5± 2.5 


10. 5± 3.0 


2.6± 1.5 







1 In units of 10~ 19 W m -2 . 



Table 5 

Observed wavelengths of the [NeV], [NeIII], and [OIII] lines in units of A 



# 


Object 


A [NeV] 3346 


A[NoV]3426 


A[NoIII]3869 


A[NoIII]3968 


A [OIII] 4959 


A[OIII]5007 


A[OI]6300 


A[OI]6363 


A[SII]6716 


A[SII]6734 


1 


PG 0026 


3833. 8±0. 4 


3923. 5±0. 4 


4431. 4±0. 4 


4545. 6±0. 3 


5679. 7±0. 3 


5734. 7±0. 3 










2 


PCj 0052 


38b0.4±0.6 


one A T_l_n r 

3y54.7±0.5 


4465. 9±0. 4 


4579.9±0.o 


5724. o±0. 3 


5779. y±0. 1 


72b9.0±2.0 


7350. 4±1. 5 


77o3.7±O.D 


77b8.D±0.5 


3 


RXJ0334 


3460. 9±0. 4 


3544. 1±0. 3 


4003. 5±0. 4 


4105. 9±0. 6 


5131. 4±0. 3 


5180. 9±0. 3 


6517.9±0.4 


6583. 4±0. 5 


6949. 3±0. 4 


6964. 0±0. 4 


4 


RXJ0337 


4284. 7±0. 5 


4387.5±1.0 


4954. 6±0. 4 


5082. Oil. 


6351. 3±0. 8 


6412. 5±1. 5 


8071.7±1.0 








5 


RX.I0354 


3461. 7±0. 6 


3545. 9±0. 6 


4005. 8±0. 5 


4108. 1±0. 5 


5134. 2±0. 2 


5183. 9±0. 3 


6522. 1±0. 4 


6586. 9±0. 4 


6952. 7±0. 0.5 


6967.7±0.5 


6 


RXJ0751 


3673.2±1.0 


3762.8±1.0 


4252. 0±0. 5 


4361. 1±0. 4 


5450. 3±0. 6 


5503. 2±0. 7 


6925. 6±0. 5 


6993.8±2.0 


7382.0±2.0 


7395.1±2.0 


7 


MS 0754 


3665. 7±0. 4 


3753. 2±0. 3 


4236. 8±0. 5 


4345. 6±0. 5 


5431. 6±0. 5 


5484. 2±0. 3 






7356. 9±0. 7 


7372. 8±0. 8 


8 


RXJ0836 


4194.8±0.6 


4295. 9±0. 5 


4852. 1±0. 5 


4976. 8±0. 5 


6219. 8±0. 4 


6279. 9±0. 5 


7904. 8±1. 5 


7989. 2±1. 5 


8424.8±1.0 


8443. 5±0. 8 


9 


MS 2128 


3631.8±1.0 


3720. 6±0. 7 


4201. 4±0. 6 


4309. 4±0. 9 


5385. 4±0. 4 


5437.5±0.5 


6842. 4±0. 4 


6911. 3±0. 4 


7294. 2±0. 5 


7309. 1±0. 8 


10 


PKS2135 


4015.9±1.0 


4112. 7±. 3 


4644. 5±0. 4 


4763. 4±0. 5 


5953. 0±0. 4 


6010. 6±0. 4 


7363.6±1.0 








11 


RXJ2256 


3563.7±1.0 


3649. 2±0. 4 


4121. 6±0. 4 


4227. 3±0. 5 


5282. 8±0. 6 


5333. 8±0. 5 


6712. 0±0. 7 


6779.9±1.0 


7154. 8±0. 4 


7169. 8±0. 4 


12 


PG 2304 




3572.2±1.0 


4033. 9±0. 5 




5170. 5±0. 3 


5220. 5±0. 3 


6567.5±0.6 


6633. 8±0. 7 


7001. 0±0. 9 


7017.3±0.8 


13 


RXJ2325 


3748. 6±0. 7 


3837. 1±0. 4 


4334. 2±0. 4 


4445. 8±0. 5 


5555. 5±0. 4 


5609. 3±0. 4 


7058. 7±0. 7 


7130.9±2.0 






14 


PG 2349 






4541. 8±1. 3 


4660.5±1.1 


5822. 3±0. 4 


5878. 4±0. 3 


7394. 6±1. 5 


7470. 1±1. 5 


7884. 9±0. 8 


7903. 5±0. 8 



Table 6 

Redshift-ordered Ratios of a 2 (i)/a 2 (0) measured from the [NeV], [NeIII], [OIII], [OI], and 

[SII] LINES, AND THE WEIGHTED AVERAGE OF ALL LINE PAIRS. 



# 


Object 




z 




a 2 (t) 


([NeV]) 






;[Noiii]) 




a 2 (t) 


([OIII]) 




a 2 (t; 

Q 2 (0 


-([OI]) 


2,., 

fl^([SH]) 


Average 




2 (t) 
2 (0) 


3 


RXJ0334 





.03478 


1 


.0054 


± 


0.0069 


1 


.0024 


± 


0.0071 





.9980 


± 


0.0086 





.9975 


± 


0.0097 


0.9883 ± 0.0380 


1 


.0007 


± 





.0040 


5 


RXJ0354 





.03536 


1 


.0171 


± 


0.0085 


1 


.0009 


± 


0.0069 


1 


.0015 


± 


0.0073 





.9863 


± 


0.0087 


1.0079 ± 0.0475 


1 


.0018 


± 





.0042 


12 


PG 2304 





.04265 


















1 


.0005 


± 


0.0085 


1 


.0020 


± 


0.0139 


1.0876 ± 0.0755 


1 


.0067 


± 





.0085 


11 


RXJ2256 





.06529 


1 


.0034 


± 


0.0126 


1 


.0051 


± 


0.0061 





.9988 


± 


0.0153 


1 


.0041 


± 


0.0181 


0.953 ± 0.0369 





.9999 


± 





.0057 


9 


MS 2128 





.08600 


1 


.0223 


± 


0.0141 


1 


.0074 


± 


0.0101 


1 


.0009 


± 


0.0123 


1 


.0000 


± 


0.0082 


0.9544 ± 0.0604 


1 


.0042 


± 





.0057 


7 


MS 0754 





.09533 





.9983 


± 


0.0057 


1 


.0064 


± 


0.0065 


1 


.0019 


± 


0.0111 










1.0097 ± 0.0675 


1 


.0023 


± 





.0046 


6 


RXJ0751 





.09914 


1 


.0199 


± 


0.0161 


1 


.0056 


± 


0.0059 


1 


.0042 


± 


0.0175 





.9776 


± 


0.0302 


0.8292 ± 0.1790 


1 


.0023 


± 





.0068 


13 


RXJ2325 





12033 





.9875 


± 


0.0090 


1 


.0009 


± 


0.0057 


1 


.0019 


± 


0.0105 


1 


.0152 


± 


0.0298 







.9987 


± 





.0048 


1 


PG 0026 





14537 





.9788 


± 


0.0062 


1 


.0099 


± 


0.0044 


1 


.0019 


± 


0.0077 















.9982 


± 





.0033 


2 


PG 0052 





15439 


1 


.0214 


± 


0.0095 


1 


.0005 


± 


0.0056 





.9945 


± 


0.0057 


1 


1041 


± 


0.0341 


0.897 ± 0.0426 


1 


.0040 


± 





.0044 


14 


PG 2349 





17404 










1 


.0240 


± 


0.0147 





.9969 


± 


0.0089 


1 


.0134 


± 


0.0285 


1.1020 ± 0.0713 


1 


.0139 


± 





.0087 


10 


PKS2135 





.20048 


1 


.0080 


± 


0.0109 


1 


.0033 


± 


0.0054 


1 


.0011 


± 


0.0983 












1 


.0038 


± 





.0046 


8 


RXJ0836 
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Fig. 1— Schematic diagram of the [OIII], [Nelll], and [NeV] (left panel), and the [Oil] and [SII] (right 
panel) line doublets 
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2.— Optical spectra of the [NeV], [Nelll], and [OIII] regions 
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3. — Optical spectra of the [01] and [SII] line regions 
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Fig. 4.— Redshift vs. a 2 {t)/a 2 {0) 



20 



